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ABSTRACT 

This paper presents Australia Telescope Compact Array Hi synthesis observations of 
the weak blue compact dwarf (BCD) galaxy NGC 2915. It is shown that NGC 2915 has 
the Hi properties of a late type spiral galaxy (Sd - Sm), including a double horn global 
profile, and Hi spiral arms. The Hi extends out to over five times the Holmberg radius, 
and 22 times the exponential scale length in the B band. The optical counterpart 
corresponds to a central Hi bar. The Hi distribution and kinematics are discussed 
in detail. A rotation curve is derived and fitted with a mass model consisting of a 
stellar disk, a neutral gas disk, and a dark matter (DM) halo. The DM halo dominates 
at nearly all radii. The total mass to blue light ratio, M. T jL^ = 76 within the last 
measured point. Thus NGC 2915 is one of the darkest disk galaxies known. The complex 
Hi dynamics of the central region results in a high uncertainty of many of the fitted 



parameters. Nevertheless it is clear that the core of the DM halo is unusually dense 
(p ps 0.1.Mopc~ 3 ) and compact (R c ~ 1 kpc). The neutral gas component, with mass 
M g = 1.27 x 10 9 A4 & is probably more massive than the stellar disk. Split and broad H I 
lines (velocity dispersion ps 35kms _1 ) are seen in the central region. Pressure support 
is probably significant, and it is not clear whether the core is in equilibrium. Beyond the 
optical disk the average Hi line of sight velocity dispersion is 8 kms -1 , which is normal 
for disk galaxies. NGC 2915 does not obey the Tully-Fisher (1977) relation, being 
underluminous for its V^ot — 88 kms -1 by a factor of nine. It also does not obey the star 
formation threshold model of Kennicutt (1989), when only the neutral gas is considered. 
A simple Hi surface density threshold of £m,crit ~ 10 21 cm~ 2 adequately describes 
the location of current star formation. Although the Hi properties of NGC2915 are 
extreme relative to normal galaxies they appear less extreme in comparison to other 
BCDS, which have similar radial profiles of Hi density and velocity dispersion, and Hi 
extending well beyond the optical disk. 
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1. Introduction 

The quest for understanding dark matter (DM) 
halos of galaxies using Hi probe is pre- 

dominantly weighted towards spiral and dwarf 
irregular galaxies. The Hi rich blue compact 
dwarf (BCD; Thuan and Martin, 1981) class has 
been largely ignored, and is under represented 
in Hi synthesis studies. This may be because 
their small optical angular sizes, often smaller 
than typical synthesized beam sizes, suggests 
that they will not be well resolved at A21 cm. 
Also, only 8% of the BCDs in Thuan and Mar- 
tin's sample have double horn profiles indicative 
of extended rotating disks, suggesting that they 
are not ideal candidates for DM studies. Some 
important Hi synthesis observations of BCDs, 
and similar galaxies, have been made (e.g. Taylor 
et al, 1993; 1995; 1996; Hunter et al, 1994) but 
their primary aims were other than to search for 
DM. 

Here we present Hi observations of NGC 2915 
whose optical properties were presented by Meu- 
rer et al. (1994; paper I). There it was shown to 
be a weak BCD, that is its integrated star for- 
mation rate is only 0.05 Mq yr^ 1 . Most of this 
star formation is near the very center of the gal- 
axy with some enhanced star formation along the 
major axis to the SE of the center. NGC 2915 re- 
solves into stars; the brightest of these yields the 
distance D = 5.3 ± 1.3 Mpc. At this distance 1" 
corresponds to 26 pc, and 1' corresponds to 1.54 
kpc. Previous single dish Hi observations indi- 
cate that it has a very extended H I disk (Becker 
et al. 1988). The Hi observations presented here 
were made primarily to determine NGC 2915's 
Hi rotation curve, and from it, constrain the 
structure of its DM halo. In addition we wish 
to examine how the Hi properties of this BCD 
differ from other types of galaxies, and whether 
they give any indication of what regulates star 
formation in BCDs. 

The observations are presented in §Ej Sec- 
tion §[I] presents an overview of the Hi proper- 



ties. The dynamics of the ISM are analyzed in 
§E]. The rotation curve is measured, pressure 
support corrections applied and mass models fit- 
ted. Section SpD presents a discussion of the re- 
sults including a comparison with other galax- 
ies, a discussion of the Tully-Fisher (1977) rela- 
tion and the star-formation threshold law. Our 
conclusions are summarized in §|6]. In the Ap- 
pendix we discuss a possible interaction partner 
of NGC 2915. 

2. Observations and Data Reduction 

NGC 2915 was observed with three configura- 
tions of the Australia Telescope Compact Array 
using all six antennas. However, baselines with 
antenna 6, (which gives the longest baselines - up 
to 6 km) were discarded since they produced in- 
significant correlation amplitudes. The dual po- 
larization AT receiver was employed with the cor- 
relator set to 256 channels per polarization, and 
with each channel separated by 15.7 KHz (3.31 
km/s). The observing logs for the three config- 
urations are given in Table [l]. The integration 
time was 20 s per visibility measurement. Some 
data in each run were lost to equipment failure 
and weather. 

The data were reduced using standard soft- 
ware in the AIPS package (ATNF version). The 
relevant tasks used are noted here in parenthesis. 
Bad or suspect data were edited out (TVFLAG, 
SPFLAG). Temporal gain and phase drifts were 
calibrated (CALIB) using the secondary calibra- 
tor, 0906-682 which was observed at 50 min. in- 
tervals for 3-6 min. Spectral variations in the 
calibration were determined from the observa- 
tions of 1934-638 and 0407-658 (BPASS). The 
absolute flux level was set by the primary cali- 
brator 1934-638 (SETJY, GETJY; flux reference: 
Walsh, 1992). The calibration was checked for 
the 1.5 km array data set using 0407-658 which is 
also a flux standard (Walsh, 1992). The calibra- 
tion agreed to within 2% with the adopted 1934- 
638 calibration. The calibrations were applied 
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and the two polarization channels combined to 
form Stokes "I" data sets (SPLIT). The contin- 
uum was fitted and subtracted in the UV plane 
(UVLSF) and the resultant continuum and line 
data sets were averaged in time (UVAVG) to one 
minute per visibility. The line data sets were 
shifted to the heliocentric rest frame (CVEL), 
and the data from the three runs were com- 
bined (DBCON). The data were imaged and 
"cleaned" (MX; Schwab, 1984; see also Clark, 
1980; Hogbom, 1974) to about the noise level 
per resultant channel to make both uniform and 
natural weighted data cubes (hereafter UN and 
NA respectively). For the NA cube, two spectral 
channels at a time were averaged at the imaging 
stage. The UN data was imaged at full spectral 
resolution. The beam sizes at the 50% level are 
W50 = 45" (circular) for the NA data cube and 
W 50 = 27" x 23" (PA = 0°) for the UN cube. The 
final cubes were made by reconvolving the clean 
component images with two-dimensional Gaus- 
sian beams having the above sizes. Planes of the 
resultant NA data cube are shown in Fig. [l|. 

From each cube, maps of total intensity, mean 
velocity, and line broadening were constructed 
from the zeroth, first, and second moments of the 
data cubes (w.r.t. velocity; MOMNT). The mo- 
ment maps and data cubes were then corrected 
for the fall-off in primary beam response with 
distance from the pointing center (PBCOR). Ta- 
ble H summarizes some of the properties of the 
NA and UN data sets including the beam size, 
pixel size and noise level in the resultant data 
cubes. 

In order to check the moment analysis results, 
single Gaussian fits to the UN data cube were 
made (XGAUS). As a further check, and to ex- 
amine the line profiles in detail the UN data cube 
was block averaged into 30" x 30" pixels, and all 
the resulting profiles were examined and fit with 
a multiple Gaussian fitting algorithm using the 
IRAF program SPLOT. 



3. Overview of HI properties 

A cursory examination of the Hi properties 
of NGC 2915 shows that they are very differ- 
ent from its optical properties. Figure I (plate 
XXXX) shows three Hi surface density maps, 
plotted as grey scales. Panel a shows the NA mo- 
ment map (45" beam); panel b, the UN moment 
map (27" x 23" beam), and panel c, the XGAUS 
peak amplitude map. NGC 2915 clearly has a 
spiral morphology, and this is most readily ap- 
parent in the NA map. There appears to be two 
arms which can be traced back to the ends of cen- 
tral bar structure. Figure || shows a contour plot 
of the NA moment map, while the UN moment 
contours are shown in Fig. ||| (Plate XXXX) 
overlaid on an / band image obtained with the 
Anglo-Australian Telescope. This clearly shows 
that the Hi, including the spiral arms, extends 
well beyond the optical extent of NGC 2915. 

The prominent central bar is about 2.5' (3.9 
kpc) long. It is most clearly seen in panel a 
of Fig. ^. The velocity field, analyzed in §4.1., 
shows the signature of an oval distortion demon- 
strating that this is a real bar, and not a highly 
inclined disk. The bar is resolved into two clouds 
in the UN data as can be seen in the bottom 
two panels of Fig. ^ and Fig. 4[ Each of these 
clouds has a neutral gas mass^ of ~ 2.5 x 10 7 
Mq (1.4 x 10 7 Mq if a local background is sub- 
tracted) which is the majority of the mass in the 
bar. They are located 17" NW and 36" SE of 
cluster 1 (paper I). They encompass many of 
the embedded objects discussed in paper I and 
bracket the brighter central clusters as well as 
some of the peculiar structure on the SE side of 
the galaxy. The faint "jets", first mentioned by 
Sersic et al. (1977) (cf. Fig. 2 paper I) project 
outwards from the core of these clouds. The 
XGAUS map in panel c of Fig. |2| is very noisy, 
and clearly shows only the highest S/N features. 
In this map the morphology is reminiscent of a 



1 The neutral gas mass is the HI mass multiplied by 1.33 
to correct for the neutral He content. 
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spiral with a bar interior to a detached ring. 

The optically visible portion of NGC 2915 co- 
incides in position, orientation, and size to the 
central bar. Thus while the exponential optical 
surface brightness profile (for R > 35") implies 
a disk structure, the Hi data suggests that the 
stars are in a bar configuration. This interpre- 
tation is consistent with the optical morphology 
(paper I) which shows enhanced star formation 
primarily in the center of the galaxy, but also 
along the optical major axis, or the bar ridge line. 
Such a distribution of star formation is common 
in barred galaxies (Priedli & Benz, 1995, and ref- 
erences therein). 

There is a large "hole" ~ 130" x 180" (3.4 x 4.6 
kpc) located 310" (8.9 kpc) to the east of the gal- 
axy's center. Examination of a position-velocity 
cut through this structure does not reveal split 
line profiles, indicating that it is probably not a 
kinematic shell, but an interarm clearing. The 
Hi in NGC 2915 does not appear to be very 
porous, as is found to be the case in many nearby 
disk galaxies, such as HoII (Puche et al., 1992). 
This may be largely a result of resolution; the 
HoII data have a resolution of 66 pc, and typical 
bubble sizes are a few hundred pc, while the lin- 
ear resolution of our data is ten times worse (640 
pc). 

The velocity field from the NA and UN mo- 
ment 1 maps are shown in Figs. [| and |6] (Plate 
XXXX) respectively. They clearly show the ex- 
pected pattern of a rotating disk. The disk must 
be warped as indicated by the twisting of the 
kinematic major and minor axes. Thus although 
there are closed velocity contours, it is not read- 
ily apparent whether these are due to a declining 
rotation curve, or the warped disk. The rotation 
curve is derived in detail in §4.1.. The NA and 



UN global velocity profiles are shown in Fig. [?]. 
These were made from the data cubes by sum- 
ming the flux in each channel within a "blanking 
aperture". The aperture was defined by eye us- 
ing the total intensity map as a guide. The twin 



horned nature of the profile is the classic profile 
of a rotating disk. The bump near V sys is due to 
the large central concentration of Hi. 

Figure |8| shows the line of sight velocity dis- 
persion, b, map from the XGAUS profile fitting. 
It shows that b increases by several tens of km s^ 1 
towards the center of NGC 2915. This is a much 
larger gradient than typically observed in spi- 
ral and irregular galaxies (Kamphuis, 1992; Lo 
et al, 1993). These maps indicate that pressure 
support may be important to the dynamics of 
NGC 2915. The line profiles are discussed in de- 



tail in § [4. 2. | , while the pressure support, or asym- 
metric drift, corrections to the rotation curve are 
discussed in 



4.3. 



Table Q presents some properties measured 
from our data. From the global Hi profiles we 
measured the total line flux / SdV, the corre- 
sponding Hi mass, Mm, the velocity widths at 
50% and 20% of the maximum intensity, W50 
and W20 respectively, and the systemic veloc- 
ity V sys . That J SdV from the UN and NA 
global profiles agree so well, indicates that the 
UN data is not missing much short spacing flux 
compared to the NA data. The Hi mass, A4hi = 
9.6 ± 0.2 x 10 8 M Q from the NA data, is signifi- 
cantly higher than that derived by Becker et al. 
(1988), Mm = 7.0 x 10 8 (after adjusting to our 
adopted distance) using a low resolution H I map 
derived from Parkes 64m observations. We sus- 
pect that the flux difference is due to the lower 
signal to noise ratio of their observations and per- 
haps inadequacies in their beam profile correc- 
tion. The systemic velocities measured from the 
global profile, V sys (global), and the tilted ring 
analysis, V sys (dynamical) (see §4.1,), agree well 



with each other, and also with the optical ve- 
locity V sys (optical) = 468 db 5kms" 1 reported in 
paper I. 

From the total intensity maps we measure 
Rm, the radius where the face-on Hi column 
density N(HI) = 5 x 10 19 cm -2 , and Ani(max), 
the maximum Hi face-on column density. The 
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correction to face-on is for the mean inclination 
i = 59°, as derived below. The remaining prop- 
erties in Table ^ compare the Hi and optical 
properties (paper I). .Mm is compared to the B 
band luminosity, L b , and Rm is compared to the 
Holmberg radius i?Ho and B band exponential 
scale length a^ 1 . The resultant ratios quantify 
what is already clear: NGC 2915 is gas-rich and 
very extended in Hi. 

4. ISM Dynamics 

4.1. Tilted ring analysis 

The Hi rotation curve was derived using the 
now standard tilted ring algorithm, ROTCUR 
(Begeman 1989), in a manner similar to that out- 
lined by Martimbeau et al. (1994). The UN and 
NA first moment maps, and the UN-XGAUS ve- 
locity field were analyzed separately, using ring 
widths of 25", and 45" for the UN and NA data 
respectively. There are up to six parameters for 
each ring; the central coordinates X c , Y c , the sys- 
temic velocity V sys , the inclination i, the position 
angle (ft of the kinematic major axis (receding 
side) , and the rotation velocity, V ro t ■ The param- 
eters were derived using five iterations of the al- 
gorithm. In the first iteration all parameters are 
allowed to be free. This iteration gives a feel for 
how the parameters vary with radius R. In the 
second iteration, i and (ft are held fixed while the 
other parameters are allowed to vary. The aim of 
this iteration is to find the best average value of 
X C ,Y C , and V sys , which are then held fixed at a 
constant value for all subsequent iterations. The 
resultant mean V sys values are reported in Ta- 
ble [2| as V sys (dynamical) . In the third iteration, 
new radial profiles of i, (ft, and U rot are derived. 
Smooth curves are drawn through these profiles 
and the interpolated values of these parameters 
are used as the initial guess for these parameters 
in the fourth iteration, which is meant to test 
how stable the solutions are to small changes in 
the initial guesses. 



It was found that the NA-MOMNT and UN- 
XGAUS solutions are stable to small perturba- 
tions. However the inner three rings of the UN- 
MOMNT field (R < 90") are not. No unique 
solutions could be found. This is due to the 
kinks in the central isovelocity contours seen in 
Fig. p. Therefore we adopt the UN-XGAUS re- 
sults for R < 90", the UN-MOMNT results for 
90" < R < 390" and the NA results for R > 390". 
Finally each half (approaching and receding) of 
the galaxy was analyzed separately. This final 
step allows the uncertainty of the parameters to 
be estimated from the level of asymmetry in the 
velocity field; the errors in the parameters are 
taken to be the absolute difference between the 
value for the full ring and either the receding or 
approaching half, whichever is larger. The V ro t 
errors estimated by this method are always larger 
than or equal to the formal errors from the least 
squares fitting. The formal errors in i and (ft are 
occasionally larger than the asymmetry of field 
errors. In these cases the former is adopted. 

The results in terms of i, (ft and V ro t for each 
data set are shown in Fig. ||[ The adopted fit pa- 
rameters are tabulated in Table ||. The Hi disk 
shows a strong warp in both (ft {A(ft = 40°) and 
i (Ai = 25°). The warp is strongest, especially 
in i, for R < 150", with relatively mild warp- 
ing (Acft = 15°, Ai = 4°) beyond this radius. 
The adopted rotation curve shows a rapid rise 
to y rot ~ 80 kms" 1 at R = 202.5". It then re- 
mains approximately constant until rising again 



for R > 290". Figure 10 shows a position velocity 
cut through the NA data cube along the major 
axis {(ft = —70°). The second rise can be seen di- 
rectly on both sides of the galaxy, demonstrating 
that it is real. 



Figure 11 shows the residual maps from the 
ring fits. Local deviations from the fit have 
\AV r \ ~ 15 kms -1 . The residuals in the in- 
ner portion show rings of mostly negative and 
mostly positive residuals at R ~ 65" and 150" 
respectively. This indicates that V sys varies with 
R, suggesting that the inner portions of the gal- 
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axy are not in equilibrium. Kamphuis (1992) 
finds V sys variations in normal spiral galaxies, 
but at large radii. In the center where V TOt is 
low these deviations are large enough to cause 
the isovelocity contour twists and the problems 
in the UN-MOMNT fit noted above. The resid- 
uals along the minor axis have opposing signs up 
to ~ ilOkms -1 on either side of the center (+ to 
the NE, — to the SW). This is a signature of an 
oval distortion^: the kinematic major and minor 
axes are non-orthogonal (Bosma, 1981). The ef- 
fect can also be discerned in Figs HH The distor- 
tion probably arises from the central bar, and is 
strongest for the rings with 200" ~ R ~ 450". 
It shows that the orbits are not strictly circular 
for these rings, although one should bear in mind 
that the residuals are small compared to the am- 
plitude of VJ-ot ~ 80 kms" 1 at these radii. 

4.2. Line profiles 

The results of the automatic Gaussian fitting 
are shown in Fig. ^ which gives a map of b values 
from XGAUS. The b map is not symmetric about 
the center of NGC 2915. There is stronger line 
broadening on the SE side than the NW. The 
broadening is centered on the SE H I cloud (see 
Fig. |). 

XGAUS was set to make only single Gaussian 
fits. The SPLOT measurements were designed 
to test the XGAUS results and measure complex 
line profiles. The position of the splot fits are 
shown in Fig. O, with the positions of profiles 
best fit by multiple Gaussians highlighted. The 
H I profiles are predominantly single peaked. The 
main region of split lines is confined to within 2' 
of the center, which is also where the lines be- 
come the broadest. Profiles at the position of 
both Hi clouds are split. Thus the total Hi 
width distribution (including splitting) is a bit 
more symmetric than the b map shown in Fig. pi. 



The average line split is AV r = 41 kms . A few 
representative profiles and their SPLOT fits are 



shown in Fig. 13 



Could the broad profiles be due to beam smear- 
ing of the observed steep rotation curve? This 
is not likely, because the knee in the rotation 
curve is four to six beam widths from the cen- 
ter, and thus the central velocity gradient is not 
that steep. Beam smearing only shows its ef- 
fects over one to two beam widths. We verified 
this by simulations of the data including the ef- 
fects of beam smearing on the UN data cube. A 
change in VJ-ot of 80 kms -1 in one beam width 
can result in b = 27 kms -1 . But we can rule 
out such a steep rotation curve, because it would 
require the knee in the rotation curve to be at 
R = 50". Only higher resolution observations 
can determine if there are large velocity gradients 
within any given beam. However any such gradi- 
ents are likely to be turbulent in nature since the 
broad (b > 20 kms -1 ) line region is eight beam 
diameters wide. 



Figure |14j compares the SPLOT b measure- 
ments with the azimuthally averaged b from the 
MOMNT and XGAUS results. It shows that the 
SPLOT and XGAUS results are in good agree- 
ment, indicating that the additional beam smear- 
ing in the 30" pixel data has not substantially 
increased b with respect to the full resolution 
UN data. The MOMNT analysis, on the other 
hand, yields b substantially lower than the Gaus- 
sian fits. This is due to an artificial bias in 
the MOMNT analysis technique; MOMNT only 
uses the data where the signal is above a cer- 
tain threshold in a smoothed version of the data 
cube. This systematically excludes the profile 
wings, and thus the results are biased towards 
lower b. 



2 Since the line widths along the minor axis are not elevated 
relative to the major axis, it is not likely that residuals are 
due to a galactic wind. 



The dotted line in Fig. 14 shows the two pixel 
resolution of the UN data cube (6 = 2.8 kms -1 ). 
Most profiles are well sampled by our data, those 
having b ~ 2.8 kms -1 mostly have low signal 
to noise ratios. A notable exception is shown 
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in Fig. 13 f. At large R, the mean b is around 
8 kms -1 , which is typical for disk galaxies (Kam- 
phuis, 1992). Figure [l^ shows that at this width 
the profiles remain well sampled. 

4.3. Asymmetric drift correction 

The broadening of the Hi profiles in the cen- 
tral region is significant compared to the mea- 
sured rotation, as can be seen by comparing 
Fig. 14 with Fig. |9[ The random motions of the 
gas may then provide significant dynamical sup- 
port in the central region. It is not clear that the 
gas there is in dynamical equilibrium. Indeed the 
asymmetry of the b map suggests that it is not. 
Nevertheless, here we will assume that gas is in 
equilibrium and account for the pressure support 
by estimating an "asymmetric drift" correction 
following Oort (1965; see his eq. 10). We assume 
that the gas velocity ellipsoid (i.e. pressure) is 
isotropic. This greatly simplifies the correction^. 
In a cylindrical coordinate system, the radial (R) 
acceleration Kr about the rotation axis z is given 
by 

>AHpb 2 }) 



-K R 



v 2 

"'rot 



R dR ' 

where p is the density. If there were no pressure 
support the orbits would be circular with veloci- 
ties V c given by 



V 2 

y c 

R 



-K R . 



The density can be written as p = H g /(2h z ), 
where h z is the vertical scale height and T, g the 
gas surface density (here taken as H I plus neutral 
He). The circular velocity is then 



Kot 2 + O'D 



(1) 



3 Binney and Tremaine (1978) treat the case of the asym- 
metric drift in collisionless systems (stars only) which have 
non-isotropic velocity ellipsoids. 



where ob is the asymmetric drift correction given 
by 



-I 



-Rb 2 



gln(E 9 ) | 2 dln(b) 



dR 



dR 



dhx{h z 
dR 



(2) 

Usually the contributions of each of the gradients 
are zero or negative, making the correction posi- 
tive. We do not have a direct measurement of h z 
but will assume dln(h z )/dR = 0. The other two 
terms can be estimated from the radial profiles 
of ln(6) and ln(E 9 ) which are shown in Fig. [l^. 
Both profiles have a core-halo structure and are 
well fit by a function consisting of a Gaussian 
centered at R = on top of a polynomial (in R) 
background: 



ln(6) 
ln(£ g ) 



2.40 + 6.57 x 10~ 4 i? 
+1.23exp(-{ J R/85.1} 2 /2) (3) 
0.812 + 1.057 x 10~ 3 i? 
-7.13 x W- 6 R 2 
+1.51exp(-{ J R/56.2} 2 /2) (4) 



where R is in arcsec, b is in kms -1 , and S 9 is the 
gas surface density in M. Q pc -2 . The fits were 
unweighted and have an rms of 0.05 and 0.07 in 
ln(6) and ln(S 5 ) respectively. These fits (con- 
verted to linear units) are shown in Fig. |1| as 
well as the ctd curve resulting from differentiating 
the fits. The last two columns of Table ^ tabulate 
ob and the corrected circular velocities, V c . As 
expected the ar> corrections are quite significant 
for R < 152". We can now proceed to model the 
mass distribution, with the caveat that the Hi 
may not be equilibrium at these radii. 

4.4. Mass models 

Three component mass models were fit to the 
rotation curve. The components are the stellar 
disk, a neutral ISM (Hi) disk, and a dark mat- 
ter (DM) halo. The stellar mass distribution is 
determined from the surface brightness profiles 
in paper I, and has one free parameter M/Lb- 
The mass distribution of the Hi disk (ncluding 
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neutral He) is determined from the observations 
presented here, and there are no free parameters. 
The DM halo is assumed to have a spherically 
symmetric density distribution given by 

9 ~ 1 + (R/Rcp { ' 

where the free parameters are the central density 
po and the core radius R c . Here we adopt 7 = 2. 
For this density distribution the rotational veloc- 
ity at large R, Vqo, and halo velocity dispersion 
do are given by (Lake et al. 1990): 

= 4ttG Po R 2 c = 4.9a . 

Fitting was done with a x 2 minimization tech- 
nique. Thus the fitted parameters depend not 
only on the input data, but also the associated 
errors. 

We experimented with numerous schemes for 
assigning errors to the V c values in Table |||. 
These included adopting the V ro t errors, com- 
bining the errors with likely ao errors, and 
taking local averages of the errors over several 
rings. In the end we adopt a uniform error in 
Vc of 3 kms^ 1 . This corresponds to the average 
error in V TO t for R > 90". The other methods 
invariably assign less weight to the inner points 
than to the outer points. The x 2 minimization 
then effectively sacrifices the fit in the inner re- 
gions in order to make a marginal improvement 
at large R. Only equal weighting produces a rea- 
sonable fit at all R. So for the adopted fits we are 
no longer doing true x 2 minimization, but effec- 
tively a non-weighted least squares minimization. 
Four models are discussed. They are illustrated 
in Fig. |l6| , and their parameters are tabulated in 
Table ||[ Each panel of Fig. [l6| shows the fitted 
rotation curve and the contribution of each of the 
components. 

Model A is the best fit to the V c data. All 
three parameters are allowed to be free, with the 
only constraint that they be > 0. The resulting 
fit has M./Lb = 0.0, i.e. the disk is not needed 
to fit the rotation curve. The halo parameters 



are for a very compact and dense core, the most 
extreme of the models generated here. 

Model B is a fit to V m t instead of V c \ in effect 
it ignores the pressure support. The difference 
between Models A and B is thus indicative of 
the uncertain dynamical state of the central re- 
gion. Neglecting the asymmetric drift correction 
results in a larger, less dense core. The effect on 
the inferred core parameters is quite dramatic, a 
factor of about 1.5 in R c and more than 2 in pq. 
Since we may have underestimated ap, it is pos- 
sible that NGC 2915's halo is even more dense 
and compact than that of Model A. Like model 

A, model B yields M/L B = 0.0. 

Model C shows a fit where the DM halo is 
a non-singular isothermal sphere (e.g. Binney & 
Tremaine 1987). For this density distribution 

00 = ^RcVnGpo = -J=. 

This is the form of the DM distribution usu- 
ally favored by Carignan and collaborators (e.g. 
Puche & Carignan, 1991, and references therein). 
It is clear that model A is a better fit over all 
radii, but model C produces as good or better of a 
fit for R < 6 kpc. The difference between Models 
A and C illustrates the uncertainty due to which 
arbitrary form of the halo distribution is adopted. 
The isothermal sphere and the analytic form of 
eq. H have significantly different scaled density 
profiles over their first few R c . Thus model C 
has very different parameter values than A and 

B. It has a much less dense core, and it is the 
only one of the three with a non-zero Ai/Ls as 
the best fit. 

Models A-C give the best fit value with all 
parameters free. The resulting limits on M/Lb 
are lower than one would expect for the red 
color of the outlying diffuse stellar population 
in NGC 2915. For Model D we return to the 
analytic halo models given by eq. || and fix 
M/Lb = 1-2, which represents approximately 
the maximum the disk component can contribute 
to match the innermost V c data point. This mass 
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to light ratio is in the range expected for the 
colors of NGC 2915 (1 £ M/L B Z 3; Bruzual 
and Chariot, 1993). The resultant halo parame- 
ters turn out to be identical to model B. In this 
case the disk contribution to V C (R) is of the same 
order as the <td corrections. Likewise models of 
type C, but with no <td corrections, produce neg- 
ligible disk contributions. Thus the strength of 
the disk component critically depends on the <jd 
corrections. Fits using the R band light profile 
(paper I) are not significantly better in fitting 
the rotation curve than the B profile used here. 
Model D is our preferred fit, since the limits on 
M/Lb are implausibly low for models A and B, 
and model D is clearly a better fit than model C. 

Note that models A, B and D have very sim- 
ilar values Voo with low uncertainties. This is 
because Voo essentially determines the level of 
the flat part of the rotation curve, which is con- 
strained much better than the inner portions. 
Voo is lower for model C because the isother- 
mal sphere reaches a maximum V c of 1.1 Voo, an d 
does not approach Voo until R ~ 10i? c (Binney 
&; Tremaine 1987). 

None of the models can produce the exact 
shape of the rotation curve, especially the up- 
turn in V TO t at large R. We performed fits where 
7 in eq. || is a free parameter. These yield 
7 = 1.93 and are not a significant improvement 
over Model A. The ±8 kms -1 modulations in the 
V c profile occur where the b profile and residual 
maps show significant non-circular motions. One 
should not over-interpret them as necessarily re- 
flecting structure in the mass distribution. 

5. Discussion 

5.1. Comparison with other galaxies 

NGC 2915 is an intriguing galaxy because its 
optical properties are those of a weak BCD (pa- 
per I) while its H I morphology and global profiles 
are those of an Sd-Sm disk galaxy (cf. Roberts 
and Haynes, 1994; Shostak, 1978). It has both a 



bar and a strong warp. These are normal prop- 
erties of spiral galaxies (Holmberg 1958; Bosma 
1991). NGC 2915 appears to be a two armed sys- 
tem. The arm structure is not flocculent like the 
transient features seen in simulations of shearing 
disks (e.g. Gerola & Seiden, 1978), despite oc- 
curring in the flat (shearing) part of the rotation 
curve. The arms are therefore probably due to a 
spiral density wave. Although the kinematic sig- 
nature of such a wave is not apparent, this may 
be due to the quality of the data (low signal to 
noise, insufficient resolution). 

NGC 2915's H I extends well beyond R Uo which 
is common for galaxies with high Mm/ Lb ra- 
tio (Huchtmeier & Seiradakis, 1985; Hoffman 
et al., 1993; van Zee et al. 1995). Its X 5 pro- 
file has a core-halo structure. This is less com- 
mon for spiral galaxies which typically exhibit 
an Hi deficit or plateau in the center (Wevers 
et al. 1986; Warmels, 1988), although similar 
cores have been found in irregular galaxies (e.g. 
NGC 55; Puche et al, 1991). NGC 2915's veloc- 
ity dispersion profile also stands out for the high 
value, b ~ 40 kms" 1 , in the center. This is much 
higher than seen in normal (non-bursting) disk 
galaxies, which typically reach a central maxi- 
mum of b ~ 15 kms" 1 , dropping to 6 to 8 kms -1 
at large radii (e.g. Dickey et al, 1990; see also 
comparisons of Kamphuis, 1992; Lo et al. 1993). 
Note that although NGC 2915 has a normal 6 = 8 
kms -1 beyond the optical radius, it is unclear 
how it and other galaxies with extended H I disks 
maintain this velocity dispersion beyond any ap- 
parent heat sources. 

The Hi properties of NGC 2915 are less ab- 
normal compared to the relatively few BCD and 
amorphous galaxies imaged in Hi. These galax- 
ies typically contain extended Hi disks well be- 
yond the optical radius, core-halo Shi profiles 
and b profiles reaching 25 to 40 kms -1 in the 
center (Taylor et al. 1994; 1995; Hunter et al. 
1994; Meurer 1994; Brinks 1988; Viallefond & 
Thuan, 1983). The properties of NGC 2915 
are also reminiscent of the starburst IBm gal- 
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axy NGC 4449 which has an Hi disk over ten 
times larger than its optical size (Bajaja et al, 
1994). This is not to say that NGC 2915 is a 
typical BCD. None of the BCDs that have been 
imaged at Hi have spiral arms. Only 8% of 
BCDs have double horn profiles and the large gas 
to stars ratio M m /L B = 2.7 of NGC 2915 also 
stands out amongst BCDs, which have a median 
M m /L B = 0.2 (Thuan & Martin, 1981). 

Less is known about the DM content of BCDs. 
This paper presents the first mass model decom- 
positions of a BCD. When Hi extends well be- 
yond i?Ho and the rotation curve is flat, the result 
is a high Mt/Lb ratio. This is certainly true for 
NGC 2915 which has M T /L B = 76 (model D) 
compared to DD0154 with Mt/L b = 74 (Carig- 
nan and Beaulieu, 1989; Hoffman et al. 1993) 
and DDO170 with Mt/L b = 57 (Begeman et al. 
1991; Lake et al. 1990). Thus NGC 2915 is one 
of the darkest disk galaxies known. We note 
that the Mt/Lb values for giant ellipticals and 
dwarf spheroidals are similarly large (e.g. Grill- 
mair et a/., 1994; Mateo et al. 1993). 

It is not just the extent of NGC 2915's Hi 
that distinguishes its remarkable DM properties. 
The DM core also appears to be extraordinarily 
dense with pq « 0.1 Mq pc~ 3 which is about an 
order of magnitude denser than typically found 
in disk galaxies (Skillman et al. 1987; Puche & 
Carignan, 1991; Begeman et al. 1991), but in the 
range of that found in dwarf spheroidal galaxies 
p = 0.07 to 1.5 M Q pc~ 3 (Mateo et al, 1993). 

Neither the core size, R c ~ 1 to 3 kpc, nor 
Voo = 73 to 90 km s _1 (depending on model form) 
are particularly distinguishing; similar values of 
R c are found in dwarf irregular galaxies, while 
low to moderate mass spirals have similar val- 
ues of Voo (Puche & Carignan 1991; and refer- 
ences therein). But having both a high Voo and 
short R c is rare. One exception is NGC 5585. 
Cote et al. (1991) fit a mass model, including 
an isothermal DM halo, to its rotation curve 
and find R c = 2.8 kpc, Vqo = 76 kms^ 1 , and 



po = 0.060.M©pc 3 , i.e. nearly identical to our 
model C. But, whereas NGC 2915 is a BCD, 
NGC 5585 is a Sd IV spiral with R Uo = 6.5 kpc, 
over twice as large at optical wavelengths as 
NGC 2915. This striking difference suggests that 
the optical morphology of galaxies is independent 
of the properties of their DM halo. NGC 5907, 
with po w 0.08 - 1.9 .M© pc~ 3 (Sackett et al, 
1994) also has one of the highest inferred po f° r a 
disk galaxy. But in this case R c is not well con- 
strained by the rotation curve and po refers to 
faint luminous matter, and not strictly speaking 
to, "dark" matter. 

It should be noted that most of the studies 
quoted above adopt distances based on Hq = 
75 km s _1 Mpc -1 . Although the distance to NGC 
2915 does not explicitly depend on Hq, it is tied 
to the distance to NGC 5253 derived by Sandage 
et al. (1994) who find H = 54kms~ 1 Mpc- 1 
(paper I). We reran our preferred model, D, tak- 
ing D = 3.1 Mpc for NGC 2915, which is more 
appropriate for Hq = 75kms _1 Mpc _1 (here D 
is derived from V r using a Virgo-centric flow 
model, and the appropriate Dvirgo), and fixing 
Mj L B = 2.0. The results make NGC 2915 
appear more extreme: R c = 0.71 ± 0.08 kpc, 
and po = 0.31 ± 0.06M Q pc~ 3 (a = 30.5 ± 0.7 
kms -1 ). If NGC 2915 were this close it would 
have Mt/Lb = 130 at its last measured point. 
A comparison of optical and DM halo proper- 
ties of galaxies employing a uniform Hq would 
be highly desirable. 

5.2. Central Energetics 

The neutral ISM is very energetic in the cen- 
tral bar area of NGC 2915. The energy is dis- 
played in two forms: that in expanding struc- 
tures, and general turbulence. The former is 
characterized by split lines implying an expan- 
sion velocity of V exp = 0.5AV, where AV is the 
amount of line splitting. The expansion energy 
is 

E k (exp) = \M g AV 2 ; 
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we take the mass involved in the expansion to 
be A4 g = 10 8 Mq (the neutral gas mass within 
the central bar area) and Ay = 40 kms -1 , and 
derive -Efc(exp) = 4 x 10 53 ergs. This is an over- 
estimate since not all of M g is likely to be in 
expanding structures. For an isotropic velocity 
distribution the turbulent energy is given by 

E fc (turb) = -M g (b HI ) 2 , 

where (bni) = 35 kms -1 is the mean veloc- 
ity dispersion in the core. Thus -Efc(turb) = 
3.8 x 10 54 erg. The energy in turbulence appears 
to be significantly greater than that in expansion, 
although beam smearing may cause some con- 
fusion between bulk flows and turbulence. We 
adopt a total kinetic energy of 

E k = £ fc (exp) + E k (tmh) = 4.2 x 10 54 erg. 

What is the source for this kinetic energy? 
Figure [D] (plate XXXX) shows the b contours of 
Fig. |H overlaid on the Ha image of Marlowe et al. 
(1995). This shows the "smoking gun" that im- 
plicates hot young stellar populations with the 
energized neutral ISM. There are two Ha bub- 
bles, of a blister morphology. That is their ioniz- 
ing clusters are at the edges of the bubbles (near 
the galactic center), and the bubbles are expand- 
ing away from these clusters. While, the profiles 
are broad over the whole ionized region, notice 
that The protrusion in the b = 30 kms -1 con- 
tour about 25" to the north of the optical center 
clearly corresponds to the Ha bubble extending 
along the optical minor axis. The most ener- 
getic Hi is associated with the SE Hi cloud. It 
is at the edge of the SE Ha bubble which ap- 
pears to be expand into it. Thus the broad Hi 
profiles are associated with the turbulent ionized 
ISM and in particular the expanding Ha bubbles. 
The ultimate energy source is the stellar winds 
and SNe produced by the young central stellar 
populations. 

For a solar metallicity population with a Sal- 
peter (1955) IMF and a constant star forma- 
tion rate, the equilibrium ratio L mec h/ '^Hq = 1-8 



(Leitherer and Heckman, 1995), where -L m ech is 
the rate mechanical energy is produced by su- 
pernovae explosions and stellar winds. From 
the = 6.4 x 10 39 ergs _1 (paper I; Mar- 

lowe et al. 1994) we deduce -L m ech ~ 1.1 x 
10 40 ergs _1 . The crossing time of the core is 
Across ~ 3.1kpc/35kms _1 = 87 Myr (i.e. the 
length of the region divided by the typical ve- 
locity). During this time the mechanical en- 
ergy release of the central population is -Emech = 
-^mech Across = 3 x 10 55 erg, or seven times E k . 
Only 10% - 20% of L m ech is likely to end up as ki- 
netic energy, the rest being thermalized (Ostriker 
& McKee, 1988; Mac Low & McCray, 1988). So 
to a first approximation, the energy released by 
the young central population in NGC 2915 is suf- 
ficient to explain the observed kinetic energy in 
the neutral ISM. 

5.3. Tully-Fisher Relationship 

Figure [E| shows NGC 2915 and DD0154 in 
the absolute magnitude versus Hi line width di- 
agram. The tight correlation shown by the filled 
circles (Sandage and Tamman, 1976) and x's 
(Sandage, 1988), is the Tully-Fisher (1977) rela- 
tionship (TFR) for nearby galaxies. NGC 2915, 
like DD0154 (Carignan and Beaulieu, 1989), 
falls way off this correlation. This discrepancy is 
also apparent in other recent calibrations of the 
TFR (e.g. Pierce and Tully, 1988) when the data 
is adjusted to the same value of Hq. Carignan 
and Beaulieu (1989) note that DD0154 shows 
the TFR may break down at the low mass end. 
NGC 2915 shows that it does not always hold 
even for moderate masses. 

Milgrom and Braun (1988) note that the TFR 
should probably be considered a correlation be- 
tween luminous matter - defined as stars and all 
phases of the ISM - and line width, not just op- 
tical luminosity and line width. The TFR holds 
for most galaxies because the stellar mass dom- 
inates over the ISM. However, for NGC 2915 
and DD0154 the neutral gas component is more 
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significant than the stellar component. If their 
ISM were converted to stars with a M/Lb = 
0.3 — 0.4 (M/Lb)q (higher if there is a signifi- 
cant mass in molecular gas) they would fall on 
the TFR. Thus under this scenario, the discrep- 
ancy with the TFR may indicate that the disks of 
DD0154 and NGC 2915 are relatively unevolved 
compared to normal galaxies. 

There still remains the problem of the dis- 
tribution of luminous matter and the TFR. In 
a normal TFR galaxy, the stellar distribution, 
which dominates the luminous mass, is much 
more concentrated than the DM. The contribu- 
tions of each component to the rotation curve 
produces similar maximum values V max , but at 
different radii, such that the overall rotation 
curve is nearly flat at V TO t = Knax from about a 
disk scale length out to the last measured point. 
This is the so-called "disk-halo conspiracy". If 
the neutral ISM of the discrepant galaxies were 
transformed to stars in situ such that the resul- 
tant luminosity is on the TFR, the stellar distri- 
bution will be coextensive with the DM (see be- 
low) and have V mSuX well below that of the DM. In 
other words there would be no disk-halo conspir- 
acy in such cases. Thus NGC 2915 and DD0154 
would not look like normal TFR galaxies even if 
their ISM were converted into stars, unless their 
ISM were to become centrally concentrated be- 
fore star formation. Likewise, one could adjust 
the distance or dust extinction As (see paper I) 
such that NGC 2915 falls on the TFR, but there 
would still be no disk-halo conspiracy and the 
stellar M/Lb would be negligible. So even then 
NGC 2915 would be unlike any normal TFR gal- 
axy. 

The coincidence of the central b = 40 km s _1 = 
Co of the DM halo is perhaps a new twist on 
the disk-halo conspiracy (although one case does 
not make a conspiracy). We speculate on the 
meaning of this coincidence in §|6]. 



5.4. Star Formation Threshold 



Figure |19ja shows the surface density of the 
contributors to the mass of NGC 2915 from 
model D: the neutral gas, £„, the stars, £*, and 



DM halo, Sdm- Figure |19| b shows the ratio of 
the first two of these with Sdm- It is immedi- 
ately apparent that the stars do not trace the 
DM very well, but the neutral gas does. A rela- 
tively flat S 9 /Edm profile is commonly observed 
in disk galaxies, as first pointed out by Bosma 
(1981; more recent examples include Carignan & 
Beaulieu, 1989; Carignan et al. 1990). Rotation 
curves of some galaxies are well fit by scaling the 
T, g profile by a constant factor to obtain the DM 
contribution (Carignan & Beaulieu 1989; Broeils, 
1992). 



Figure 19 also shows the critical density for 



high mass star formation, S cr it (Kennicutt, 1989) 
and E cr it/SDM- E 9 is less than E cr i t by a factor 
of about 3 (range of 2 to 9), even in the center 
of NGC 2915 where high mass star formation is 
seen. When only the neutral ISM is considered, 
NGC 2915 does not obey Kennicutt's star for- 
mation law. The discrepancy between S 9 and 
S cr it is worst in the center of NGC 2915 because 
S cr jt oc 6, and this is where b is the highest. How- 
ever, even adopting a constant 6 = 5 kms" 1 pro- 
file does not bring £ C rit < T, g . This is contrary 
to what is found for a wide variety of galaxies in- 
cluding both low surface brightness galaxies (van 
der Hulst et al. 1993) and other BCDs (Taylor 
et al, 1994) which obey the Kennicutt law when 
only the neutral ISM is considered. Of course 
there may be no real discrepancy, since much of 
the ISM may be in molecular form. Star forma- 
tion in NGC 2915 occurs at R < 30" (paper I), 
where E crit /S g ps 9 and M g (R < 30") = 3.1 x 10 7 
M.q. Thus to satisfy Kennicutt's law there 
should be M mo i ecu i ar > 3 x 10 8 M.q. Only about 
1% of the so-called dark matter need be molecu- 
lar gas for the Kennicutt law to hold. 

The peak H I surface density occurs in the two 
central Hi clouds and is Shi = 2.4 x 10 21 cm -2 



13 



or after inclination correction Shi = 1-2 x 10 
cm -2 . Recent high mass star formation is seen 
near the positions of both these clouds. Star for- 
mation in NGC 2915 is consistent with "simple" 
Hi threshold schemes for high mass star forma- 
tion; for example Gallagher and Hunter (1984) 
propose that efficient high mass star formation 
occurs where Shi > £ffl,crit = 5 x 10 20 cm -2 , 
while Skilhnan (1987) similarly proposes that 
Sm,crit = 10 21 cm~ 2 . Thus although, the more 
complex star formation law of Kennicutt (1989) 
does not hold for NGC 2915, a simple threshold 
law does. 

Kenney et al. (1993) show how tidal shear can 
limit star formation via cloud destruction. In 
NGC 2915 and the BCDs in the sample of Taylor 
et al. (1994), the star forming region and indeed 
most of the optically visible galaxy is located in 
the rising portion of the rotation curve. This is 
where the rotation is almost solid body, and thus 
tidal shear is minimized. Similarly the galaxies 
in the Skilhnan (1987) study all have solid body 
rotation where their H II regions are located. 

6. Conclusions / Summary 

Our detailed Hi study of the weak BCD gal- 
axy NGC 2915 shows it to have a nature unex- 
pected from its optical appearance. The H I mor- 
phology and global properties (mass, line width) 
are those of a late type barred spiral galaxy. The 
presence of a bar is confirmed by the signature 
of an oval distortion in the velocity field. At our 
highest resolution the bar is resolved into two 
Hi clouds separated by 1.5 kpc. The Hi line 
of sight velocity dispersion in this central region 
is very high, b pa 40 kms^ 1 , and the Hi line is 
split by AV pa 40 kms -1 at the position of both 
Hi clouds. The energetics of the central region 
can be accounted for by the known young stellar 
populations there. The Hi extends well beyond 
the optical galaxy, out to five times the Holm- 
berg radius. At these radii, the dynamics are 
dominated by rotation, although there are indi- 



cations of an oval distortion in the velocity field, 
and some warping of the disk. The rotation curve 
is approximately flat for R > 4 kpc, but slightly 
rising at the last measured points (R ~ 10 kpc). 
For R < 3 kpc, random motions are significant 
and there are indications that the systemic ve- 
locity varies with radius. It is not clear whether 
the system is in equilibrium at these radii. 

The rotation curve can not be accounted for 
by the luminous components alone. Large amounts 
of dark matter (DM) are required, dominating 
at nearly all radii. Within the radius of the 
last measured point (R = 15.2 kpc) the total 
mass to blue light ratio is Mt/Lb = 76, mak- 
ing NGC 2915 one of the darkest disk galaxies 
known. The complex dynamics of the central re- 
gion, make interpretation of mass model fits to 
the rotation curve difficult. Despite this uncer- 
tainty, it is clear that NGC 2915's DM distri- 
bution has an unusually dense (central density 
po sa O.lA^epc -3 ) and compact core (radius 
R c = 0.8 to 2.5 kpc). The disk M/L ratio is 
not well constrained by our models, although it 
seems likely that the stellar component is less 
massive than the neutral ISM disk. 

NGC 2915 does not obey the Tully-Fisher 
(1977) relationship, being under luminous by a 
factor of ten for its line width (or rotational ve- 
locity). Nor does it obey the Kennicutt (1989) 
star formation threshold law. Unless there is 
considerable amount of molecular material, its 
gas density is always below the critical density 
for high mass star formation, even in the cen- 
ter where high mass stars are observed to have 
formed recently. However, it does obey simple 
critical density threshold models for high mass 
star formation (c.f. Gallagher & Hunter 1984; 
Skilhnan 1987). 

Finally we speculate on the nature of star for- 
mation in NGC 2915, a topic our observations 
shed new light on. In particular they reveal two 
key ingredients in its recipe for regulating star 
formation: a bar, and a deep central potential. 
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Friedli & Benz (1995) showed through numerical 
simulations that there is a complex interaction 
between star formation, the ISM, and the bar in 
isolated barred galaxies. The bar redistributes 
angular momentum resulting in an inflow along 
the bar. Some star formation along the bar ener- 
gizes the ISM (via the resultant supernovae ex- 
plosions and stellar winds) enough to inhibit ad- 
ditional star formation until the material reaches 
the center. The coincidence in velocity disper- 
sions, that of Hi measured in the center of the 
galaxy, and that inferred for the DM halo, sug- 
gests that the intensity of the central star forma- 
tion is regulated by the depth of the potential 
well (Bothun et al, 1986); if the star formation 
becomes too strong the ISM is thrown clear of 
the core and star formation is halted, at least 
temporarily. Some models indicate that such an 
outflow can contribute to the spin down of the 
disk resulting in further inflow (Corbelli & Sal- 
peter, 1988; Charlton & Salpeter, 1989) thus es- 
tablishing a "galactic fountain" . The stellar com- 
ponent of NGC 2915 and other BCD/amorphous 
galaxies is limited to within the turnover radius 
indicating that tidal shear probably strongly in- 
hibits star formation. The optical images of 
NGC 2915 has the size and orientation of the cen- 
tral bar, suggesting that the stellar counterpart 
is a bar, or at least that the bar axis is coinci- 
dentaly the optical major axis. A scenario of bar 
and DM regulated star formation for NGC 2915 
does not require interaction, which is convenient 
since NGC 2915 is in a low density environment. 
However, there is one possible interaction part- 
ner (Appendix) so an interaction model such as 
that sketched by Noguchi (1988) or Taylor et al. 
(1993) can not be ruled out. 
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A. Appendix: an interaction partner? 

Taylor et al. (1993; 1995; 1996) surveyed ap- 
parently isolated BCDs for companions with the 
rationale that interaction with an unseen com- 
panion may be responsible for the BCDs present 
burst of star formation. Such companions were 
successfully found in 12/21 galaxies in their sam- 
ple. We unsuccessfully searched our Hi data 
cube for faint companions. However the primary 
beam of the AT dishes is W50 = 33' which corre- 
sponds to a radius of 25 kpc. At a typical relative 
velocity of 300 kms -1 an interacting partner can 
transverse this distance in less than 100 Myr. 

We searched the NASA/IPAC Extragalactic 
Database (NED)^ for possible interaction part- 

4 The NASA/IPAC Extragalactic Database (NED) is op- 
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ners within 5° of NGC 2915. There are no galax- 
ies with cataloged velocities within 500 kms -1 
of NGC 2915's. However there is one low sur- 
face brightness object, SGC0938. 1-7623, which 
is separated from NGC 2915 by 42', a projected 
distance of 64 kpc, and has no measured veloc- 
ity. Is SGC0938. 1-7623 the interacting partner 
of NGC 2915? 

Unfortunately the answer is not yet clear. 
Corwin et al, (1985) note that it is not cer- 
tain whether SGC0938. 1-7623 is a low surface 
brightness galaxy or diffuse nebulosity associ- 
ated with our galaxy. Figure ^ shows an / 
band CCD image of SGC0938.1-7623 obtained 
using the Fl imaging system on the AAT. It was 
made from four 200s exposures obtained during 
non-photometric conditions. There are numer- 
ous foreground stars over the face of SGC0938.1- 
7623, but no obvious overabundance of sources 
that would suggest SGC0938. 1-7623 is resolved. 
The main body has a diameter of ~ 130". There 
are faint extensions to the SE and NW yielding a 
total diameter of ~ 270". The overall symmetry 
of the system suggests that it is a galaxy. 

An Hi spectrum in the direction of SGC- 
0938.1-7623 was obtained with the Parkes 64m 
radio telescope. The observations cover the range 
V r = -1000 to +4000 kms" 1 split over 1024 
channels and reach a noise level of 50 mJy There 
are no obvious Hi features in the spectrum ex- 
cept for strong galactic Hi and a narrow spike 
having V r = 300 kms -1 , W50 = 30 kms -1 , and 
/ Sdv = 0.72 ± 0.06 Jy kms -1 . This probably 
corresponds to the edge of the high velocity cloud 
HVC 300-24+274, which Morras (1982) shows to 
have narrow profiles and velocities varying over 
the range 240 ~ V r ~ 330 kms -1 . 

We conclude that there is no evidence that 
SGC0938. 1-7623 is interacting with NGC 2915. 
However, it can not yet be ruled out as an inter- 

erated by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the National 
Aeronautics and Space Administration. 



action partner. For an assumed 50 kms 1 line 
width the Hi spectrum should be able to de- 
tect A4 H i > 2.6 x 10 6 M Q with S/N > 5 at 
NGC 2915's distance. If SGC0938.1-7623 is at 
the distance of NGC 2915 it is gas poor. 
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Fig. [|. Contour plots of the channels with sig- 
nal in the CLEANed NA data cube. The con- 
tour levels are at 10, 25, 50, 75, and 90 percent 
of the peak surface brightness of 4.7 Jy beam -1 
(JVm = 1-69 x 10 20 cm- 2 ). The V r in kms" 1 of 
each channel is shown in the upper right of each 
panel. The beam size (W50 = 45") is shown on 
the lower left of the upper left panel, and again 
four panels beneath it. 

Fig. ^. Grey scale representations of Hi inten- 
sity maps, all to the same scale. Top panel is the 
NA moment map (45" resolution) . The middle 
panel is the UN moment map (25" resolution). 
The bottom panel shows the peak intensity (am- 
plitude) from the XGAUS single Gaussian fits to 
the UN data. North is up and east is to the left 
in all panels. The scalebar in the bottom panel 
is 2' long and 20" wide. 

Fig. ||. Contour plot of total Hi intensity from 
the NA zeroth moment map. Contour levels are 
at 2.5, 7.5, 15, 25, 50 75, and 90 percent of 
the peak Hi surface brightness of 3.1 Jy beam -1 
kms -1 (N m = 1.71 x 10 21 cm" 2 ). In this and the 
remaining maps the crosses mark the positions of 
fiducial stars and the star marks the optical cen- 
ter from paper I, and the beam size is shown in 
the lower left corner. 

Fig. |. 800s I band CCD image of NGC 2915, 
obtained with the AAT Fl system overlain with 
H I total intensity contours from the UN moment 
map. Contours are at 5, 13, 25, 50, 75, and 90 
percent of the peak surface brightness of 1.39 Jy 
beam -1 kms -1 (N m = 2.46 x 10 21 cm -2 ). The 
strong glow at lower right in this figure and Fig. || 
is due to a bright star just off the edge of the 
frame. 

Fig. |5[ Velocity field from the NA first moment 
map. The V r in kms -1 of the contours are indi- 
cated. 

Fig. |6[ Velocity field from the UN first moment 
map overlain on the / band image. The con- 
tour interval is 10 kms -1 . The contour pass- 
ing through the center has V r = 470 kms -1 ; the 



prominent closed contours to the SE and NW of 
the center have V r = 400 and 540 kms -1 , respec- 
tively. 

Fig. ^ Global H I V r profiles derived from the NA 
(solid line) and UN (dotted line) data cubes. The 
systemic velocity (first moment of the profiles) 
are indicated with arrows of the corresponding 
line type. 

Fig. H Line broadening, b, map from automatic 
single Gaussian fitting using XGAUS. The results 
are shown as both a grey scale and contours at 
b = 10, 20, 30, 40 and 50 kms -1 . Blank areas are 
regions of either low signal to noise or complex 
profiles which could not be fit with a single Gaus- 
sian. 

Fig. |H Results of the tilted ring analysis of the 
three velocity fields: NA-MOMNT, UN-MOMNT, 
and UN- XGAUS. The parameters shown are ma- 
jor axis position angle, </>; inclination i, and rota- 
tion velocity V^ t- 

Fig. 10 Contours through a position - velocity 



cut of the NA data cube. The cut is through the 
center of NGC2915 and at a position angle of 
70°, corresponding to the kinematic major axis. 
The X axis indicates location along this position 
angle; the tick marks are separated by 208". The 
contours are at -3.5, 3.5, 7, 14, 21, 28, and 35 
mJy/Beam. Note the "second rise" in U rot seen 
on both sides of the galaxy. 

Fig. [jT| Grey scale plus contour plots of the resid- 
ual maps from the tilted ring fits. The contours 
shown are for AV r = -20, -15, -10, -5, 5, 10, 15 
and 20 kms -1 . Panels a, b, and c show the NA- 
MOMNT, UN-MOMNT, and UN-XGAUS resid- 
uals respectively. Note that the panels have dif- 
ferent scales. 



Fig. 12 Position with respect to the center of 
NGC 2915 (shown as a star) of SPLOT profile 
measurements. Positions marked with small dots 
are well represented by single Gaussian fits, while 
double Gaussians are required for the positions 
marked with large dots. The crosses mark the 
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position of the fiducial stars shown in Figs. ^],||. 

Fig. [f3| Examples of SPLOT Gaussian fits (con- 
tinuous lines) to H I line profiles (histogram style 
lines) . The profiles were chosen to illustrate (a) a 
typical high signal/noise outer profile; (b) a typi- 
cal low signal/noise outer profile; (c) a broad sin- 
gle profile near the center; (d) a split profile near 
the center; (e) a split profile in the outer disk; 
(f) an atypical profile: a narrow core on top of a 
normal single profile. The positions of the pro- 
files in ±E, dbN offsets (in arcsec) from the center 
of NGC 2915 are: (a) -295,101; (b) -25,311; (c) 
5,11; (d) 35,-19; (e) -265,-79; (f) -325,41. 

Fig. [l4] Comparison of line of sight velocity dis- 
persion, b, measurements. Azimuthal averages 
of the UN-MOMNT, and UN-XGAUS measure- 
ments are shown as filled squares and triangles 
respectively. Individual SPLOT measurements 
with signal/noise > 3 are shown as x's, lower sig- 
nal/noise splot measurements are indicated with 
dots. 



ted line; neutral gas, long dashed line; DM halo, 
short dashed line. See the text for a description 
of the individual models. 



Fig. 15 The top panel shows the line of sight line 



broadening, b, as a function of radius, while the 
middle panel shows the gas surface density, S g . 
The filled circles show the data which are az- 
imuthal averages extracted from the UN-XGAUS 
b map (top) and UN zeroth moment maps (cen- 
ter) using concentric rings with <j) = —60.8, 
i = 59.6. The broken lines show the fits to these 
data as given in the text. The contribution of the 
b and S gradients to the total asymmetric drift 
correction, <td are shown in the bottom panel 
with the corresponding line style. The total T,d 
curve is shown as the solid line. 

Fig. |l^ The rotation curve with mass model fits. 
Triangles represent results from the UN-XGAUS 
data, circles from the UN-MOMNT data, and 
squares from the NA-MOMNT data. Closed 
symbols are V c and open symbols V ro t- The 
Holmberg radius is indicated by the arrow. The 
thick solid line in each panel shows the model 
fit, and the broken lines show the contribution 
of the individual components: stellar disk, dot- 



Fig. 17 Contours of line of sight velocity dis- 
persion b overlain on the Ha image of Marlowe 
et al. (1994). The b contours correspond to those 
shown in fig. |8|. The scale bar at lower right is 
30" long. 

Fig. [l^ The Tully-Fisher relationship for nearby 
disk galaxies compared with DD0154 and NGC 
2915. AV21 is the inclination corrected Hi line 
width at 20% of the peak value, and Mb is the 
absolute B band magnitude. The solid circles are 
data from Sandage &: Tamman (1976); the x's 
are from Table 2 of Sandage (1988). Galaxies 
with i < 40° have been excluded from both sam- 
ples and we have taken Mg = M pg + 0.1 where 
necessary. 

Fig. [l9| Upper panel: radial variation of logarith- 
mic surface densities. The neutral gas measure- 
ments (solid line) are azimuthal averages using 
the UN data for R < 10 kpc (390") and the NA 
data for R > 10 kpc. The projected DM den- 
sity (dashed line) is that derived from model D. 
The stellar profile (dotted line) is the B band 
surface brightness profile (paper I) normalized to 
the M/Lb ratio of model D. The critical density 
required for efficient star formation (Kennicutt, 
1989; dot dashed line) is derived from the rota- 
tion curve and the b fit of eq. ^[ Lower panel: log- 
arithm of surface densities normalized by Sdm- 
Line styles are the same as for the upper panel. 

Fig. H 800 s I band CCD image of SGC0938.1- 
7623, obtained with the AAT Fl system. North 
is up East is to the left and the image is 500" 
across in each dimension. 
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Table 1 
HI Observing Log 



config. 


UT date 


baseline 


on source 


name 




range a 


(hrs) 



0.375 3 May, 1992 31 - 459 11.34 
0.75 C 22 Jan, 1993 46 - 750 10.52 
1.5 D 12 Mar, 1993 107 - 1439 10.80 



a In meters, excluding baselines with an- 
tenna 6. 



Table 2 



Dataset properties 



Quantity- 


NA value 


UN value 


units 


beam W50 


45 x 45 


27 x 23 


arcsec x arcsec 


pixel size 


14 x 14 x 6.62 


10 x 10 x 3.31 


arcsec x arcsec x kms" 1 


noise/channel 


2.9 


3.0 


10~ 3 Jy Beam" 1 


JSdV 


145 


139 


Jy kms -1 


Mm 


9.58 


9.21 


1O 8 7W 


W 50 (global) 


151 


146 


kms -1 


W 2 o (global) 


170 


163 


kms -1 


V sys (global) 


471 


468 


kms -1 


V sys (dynamical) 


467 


469 


kms -1 


Rm 


14.9 




kpc 


./Vhi (max) 


0.88 


1.26 


10 21 cur 2 


M m /L B 


2.7 


2.6 


(M/L B ) G 


Rbi/Rho 
Rm/ a B l 


5.2 






22.6 
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Table 3 

Results of tilted ring analysis. 



R 

( ) 


4> 
( ) 


i 

( ) 


Viot 

(km s~ ) 


(km s~ ) 


v c 

(km s 


27.5 


-28 ± 12 


73 ±24 


29 ± 17 


26.5 


39 


52.5 


-28 ±11 


78.0 ±7.1 


37 ±22 


38.1 


53 


77.5 


-48.1 ±2.1 


63 ± 11 


53 ±27 


37.7 


65 


102.5 


-56.3 ±0.8 


74.8 ± 5.4 


65.5 ±2.0 


31.5 


72.7 


127.5 


-56.8 ±3.9 


73.1 ±9.7 


71.9 ±9.0 


24.6 


76.0 


152.5 


-55.7 ±4.2 


70.3 ±8.4 


76.8 ±6.7 


19.0 


79.2 


177.5 


-53.7 ±0.6 


62.2 ±2.3 


79.9 ± 1.3 


15.2 


81.4 


202.5 


-53.0 ±0.6 


56.0 ±2.7 


81.5 ±1.0 


12.6 


82.5 


227.5 


-56.3 ±0.5 


53.8 ±1.4 


82.1 ±4.6 


11.2 


82.8 


252.5 


-58.6 ± 1.7 


55.2 ± 1.4 


79.7 ±4.6 


10.5 


80.4 


277.5 


-61.4 ± 1.9 


53.0 ± 1.2 


79.7 ±3.7 


10.5 


80.4 


302.5 


-62.8 ±2.1 


52.2 ± 1.1 


80.2 ± 1.6 


10.8 


80.9 


327.5 


-64.0 ± 1.0 


53.1 ±1.5 


79.9 ±0.7 


11.3 


80.7 


352.5 


-65.0 ± 1.1 


53.8 ±3.1 


79.9 ± 1.5 


11.9 


80.8 


377.5 


-65.1 ± 1.0 


54.6 ± 2.3 


81.0 ±0.7 


12.5 


82.0 


412.5 


-66.3 ±0.4 


56.2 ±2.0 


81.2 ±3.2 


13.3 


82.3 


457.5 


-67.5 ±1.6 


58.4 ±0.7 


82.6 ±2.1 


14.3 


83.8 


502.5 


-68.5 ±2.2 


58.7 ±1.9 


84.2 ± 1.5 


15.2 


85.5 


547.5 


-68.8 ±0.9 


58.9 ±2.3 


87.0 ±0.5 


16.1 


88.5 


592.5 


-68.8 ±0.5 


58.2 ± 1.1 


92.3 ± 1.3 


16.9 


93.9 
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Table 4 
Mass models 



Model: 



A 



B 



D 



Fit parameters : 
M/Lb (solar units) 
R c (kpc) 
p (A4 pc- 3 ) 



-L 



- 1 ) 
-V 



o"o (km 
Voo (kms 
rms (kms 
Reduced \ 2 

At R Ko (2.9 kpc) : 
V c (kms" 1 ) 
M T (1O 9 A4 ) 
M.t/Lb (solar units) 

-A^dark/ Luminous 

Within last point (15.2 kpc) : 
V c (kms" 1 ) 
M T (1O 9 A4 ) 
Mt/Lb (solar units) 

-Mdark/ -^Luminous 



< 0.4 

o.78±°;^ 

0.24 ±0.06 
39.9 ±0.8 
88 ±2 
2.8 
1.01 



73 
3.6 
10.1 
27 

87 
27 
74 
28 



< 0.14 
1 94+ - 19 
0.10 ±0.02 
41.2 ±0.9 
91 ±2 
3.8 
1.87 



66 
3.0 
8.4 
22 

87 
27 
75 
28 



0.75 ± 0.40 
2.5 ±0.3 
0.07 ±0.02 
51.9 ±0.9 
73 ±1 
4.1 
2.14 

73 
3.6 
10.2 
8.3 

81 
23 
65 
19 



1.2 

1.23 ±0.15 
0.10 ±0.02 
41.0 ± 1.1 
91 ±2 
3.8 
1.80 

71 
3.4 
9.6 
5.2 



27 
76 
19 



Note. — Errors are 90% confidence limits and upper limits are 95% confidence limits. 
They were determined from Monte-Carlo simulations assuming uniform V c and V TO t 
errors of ±3 kms -1 . 
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